Floquet Topological Insulator in Semiconductor Quantum Wells 
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Topological phase transitions between a conventional insulator and a state of matter with topolog- 
ical properties have been proposed and observed in mercury telluride - cadmium telluride quantum 
wells. We show that a topological state can be induced in such a device, initially in the trivial 
phase, by irradiation with microwave frequencies, without closing the gap and crossing the phase 
transition. We show that the quasi-energy spectrum exhibits a single pair of helical edge states. The 
velocity of the edge states can be tuned by adjusting the intensity of the microwave radiation. We 
discuss the necessary experimental parameters for our proposal. This proposal provides an example 
and a proof of principle of a new non-equilibrium topological state, Floquet topological insulator, 
introduced in this paper. 



Topological phases of matter have captured our imag- 
ination over the past few years, with tantalizing proper- 
ties such as robust edge modes and exotic non-Abelian 
excitations [TJ [2], and potential applications ranging from 
semiconductor spintronics [3 to topological quantum 
computation The discovery of topological insula- 

tors in solid-state devices such as HgTe/CdTe quantum 
wells [51 [5J, and in materials such as Bi 2 Te 3 , Bi 2 Sn 3 [7J- 
12] brings us closer to employing the unique properties of 
topological phases in technological applications. 

Despite this success, however, the choice of materi- 
als that exhibit these unique topological properties re- 
mains rather scarce. In most cases we have to rely on 
serendipity in looking for topological materials in solid- 
state structures and our means to engineer Hamiltonians 
there are very limited. Therefore, to develop new meth- 
ods to achieve and control topological structures at will 
would be of great importance. 

Our work demonstrates that such new methods are 
indeed possible in non-equilibrium, where external time- 
dependent perturbations represent a rich and versatile re- 
source that can be utilized to achieve topological spectra 
in systems that are topologically trivial in equilibrium. In 
particular, we show that periodic-in-time perturbations 
may give rise to new differential operators with topolog- 
ical insulator spectra, dubbed Floquet topological insu- 
lators (FTI), that exhibit chiral edge currents in non- 
equilibrium and possess other hallmark phenomena asso- 
ciated with topological phases. These ideas, put together 
with the highly developed technology for controlling low- 
frequency electromagnetic modes, can enable devices in 
which fast switching of edge state transport is possible. 
Moreover, the spectral properties of the edge states, i.e., 
their velocity, and the bandgap of the bulk insulator, can 
be easily controlled. On a less applied perspective, the 
fast formation of the Floquet topological insulators in re- 
sponse to the external field opens a path to study quench 



dynamics of topological states in solid-state devices. 

The Floquet topological insulators discussed here 
share many features discussed in some previous works: 
The idea of achieving topological states in a periodic 
Hamiltonian was also explored from the perspective of 
quantum walks in Ref. QUI Also, a similar philosophy led 
to proposals for the realization of topological phases in 
cold-atom systems: a quantum Hall state using a stro- 
boscopic quadrupole field and a topological insulator 
state using a Ramann-scattering induced spin-orbit cou- 
pling [12l. Also, Ref. [13] proposed to use a circularly- 
polarized light field to induce a Hall current in graphenc. 
Another useful analogy for our work is the formation of 
zero-resistance state in Hall bars at low magnetic fields 
using RF radiation [HHT7j . In our case, it is not the 
resistance of the bulk that vanishes, but rather the re- 
sistance of the edges, which becomes finite and universal 
upon application of the light field. As we were finaliz- 
ing our draft, we also became aware of Ref. [T5J which 
proposed to use a periodic modulation in the form of a 
circularly polarized light to change the Chern number in 
the Haldane model [19]. 

DEFINITION OF A FLOQUET TOPOLOGICAL 
INSULATOR 

Let us first provide a general construction and def- 
inition for a Floquet topological insulator in a generic 
lattice model, and then discuss a specific realization: 
a HgTe/CdTe quantum well. The generic many-body 
Hamiltonian of interest is 

H(t)= ^^(k.ijc^k+h. c, (1) 

keBZ 

where k and c mj k are fermion creation/annihilation 
operators, k is the momentum defined in the Brillouin 
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zone, and the Latin indices, n, m = 1, 2, ... TV label some 
internal degrees of freedom (e.g., spin, sublattice, layer 
indices, etc.). The N x N k-dependent matrix H(k,t) 
is determined by lattice hoppings and/or external fields, 
which are periodic in time, %{T + 1) = T-L(t). 

First, we recall that without the time-dependence, the 
topological classification reduces to an analysis of the ma- 
trix function, H(k), and is determined by its spectrum. 
[2"TJl |2"T] . An interesting question is whether a topologi- 
cal classification is possible in non-equilibrium, i.e., when 
the single-particle Hamiltonian, H(k,t), in Eq. does 
have an explicit time-dependence and whether there are 
observable physical phenomena associated with this non- 
trivial topology. Consider the single-particle Schrodinger 
equation associated with Eq. 0: 

[If (k, t) - ildt] * k (t) = 0, with H(k, t) = H(k, t + T) 

(2) 

The Bloch-Floquet theory states that the solutions to 
Eq. ^ have the form * k (t) = S k {t)^ k (0), where the 
unitary evolution is given by the product of a periodic 
unitary part and a Floquet exponential 

S k {t) = A(t) exp [-iH F {k)t] , with V k {t) = A(t + T) 

(3) 

where Hp(k) is a self-adjoint time-independent matrix, 
associated with the Floquet operator [if (k, t) — ildt] 
acting in the space of periodic functions $(t) = $(i + T), 
where it leads to a time-independent eigenvalue prob- 
lem, [H(k,t) -Ud t ] $(M) = e(k)$(k,i). The quasi- 
energies e(k) are the eigenvalues of the matrix Hp(k) in 
Eq. ([3]), and in the cases of interest can be divided into 
separate bands. The full single-particle wave- function is 
therefore given by ^(t) = e~ i£t Q(t). Note that the quasi- 
energies are defined modulo the frequency u — 2-k/T. 

The Floquet topological insulator is defined through 
the topological properties of the time-independent Flo- 
quet operator Hp (k), in accordance with the exist- 
ing topological classification of equilibrium Hamiltonians 
[201 EE]- Most importantly, we show below that the FTI 
is not only a mathematical concept, but it has immediate 
physical consequences in the form of robust edge states 
that appear in a finite system in the non-equilibrium 
regime. The density profile of these hallmark bound- 
ary modes was found to be only weakly time-dependent 
[through the envelope function determined by the matrix 
V k (t) in Eq. [3]. Furthermore, we explicitly demonstrate 
within a specific model that not only topological insulator 
properties may survive in non-equilibrium, but they can 
be induced via a simple non-equilibrium perturbation in 
an otherwise topologically trivial system. Using the spe- 
cific model describing a HgTe/CdTe quantum wells, we 
show that a carefully chosen non-equilibrium perturba- 
tion may be utilized to turn the topological properties 
on and off. We discuss various methods to realize such a 
perturbation using experimentally accessible electromag- 
netic radiation in the microwave-THz regime. 



TOPOLOGICAL TRANSITION IN HgTe/CdTe 
HETEROSTRUCTURES 



Below we outline a proposal for the realization of a 
FTI in Zincblende structures such as the HgTe/CdTe het- 
erostructure which are in the trivial phase (d < d c in Ref. 
EJ). Consider the Hamiltonian 



= ( H(k) 







where 



/f(k) = e(k)/ + d(k)-<7, 



(4) 



(5) 



k = (k x ,k y ) is the two dimensional wave vector, and 
& = &y,<jz) are the Pauli matrices. The vector, d(k), 
is an effective spin-orbit field. Equations Q and ([5| rep- 
resent the effective Hamiltonian of HgTe/CdTe hetero- 
junctions [5]. In these systems, the upper block H(k) is 
spanned by states with mj = (1/2, 3/2), while the lower 
block is spanned by states with mj = (—1/2, —3/2). The 
lower block H*(— k) is the time reversal partner of the 
upper block, and in the following we shall focus our at- 
tention to the upper block of Eq. Q only. 

The Hamiltonian ([5]) is diagonalized via a k- 
dependent SU(2) rotation, that sets the z-axis for the 
pseudo-spin along the vector d(k). There are two double- 
degenerate bands with energies e±(k) = e(k) ± |d(k)|. 
Within each sub-block, the TKNN formula provides the 
sub-band Chern number [22], which for the Hamilto- 
nian ([5| can be expressed as an integer counting the num- 
ber of times the vector d(k) wraps around the unit sphere 
as k wraps around the entire FBZ. In integral form, it is 
given by 



C+ =±- 
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d 2 kd(k)- d kx d{k) x d ky d(k) 



(6) 



where d(k) = d(k)/ |d(k)| is a unit vector and the (±) in- 
dices label the two bands. 

This elegant mathematical construction yields also 
important physical consequences, as it is related to the 
quantized Hall conductance associated with an energy 
band, 



-a 



(7) 



Around the T point of the first Brillouin zone (FBZ) 
we can expand the vector d(k) as [5l [23] 



d(k) = (Ak x ,Ak y ,M - Bk 2 ), 



(8) 



where the parameters A < 0, B > and M depend on 
the thickness of the quantum well and on parameters of 
the materials. We can easily see that the Chern number 
implied by d(k) depends crucially on the relative sign 
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FIG. 1. Energy dispersion e(k) and pseudospin configuration 
— d(k) for the original bands of H(k) in the non topological 
phase (A'l/B < 0). The non-topological phase is characterized 
by a spin-texture which does not wrap around the unit sphere. 
Upon application of a periodic modulation of frequency uj 
bigger than the band gap, a resonance appears; the green 
circles and arrow depict the resonance condition. 



of M and B. Within the approximation of Eq. ([8]), far 
away from the T point, d(k) must point south (in the 
negative z direction). At the V point, d(k) is point- 
ing north for M > 0, but south for M < 0. For the 
simplified band structure, the Chern numbers are clearly 
C± = ± [1 + sign(M/B)] /2. For a generic band struc- 
ture corresponding to Eq. ([8| near the T-point, the same 
logic applies, and we can easily see that a change of sign 
in M induces a change of the Chern number, C, by 1. 

We now show that a similar non-trivial topological 
structure can be induced in such quantum wells, start- 
ing with the non-topological phase, via periodic modu- 
lation of the Hamiltonian, which allows transitions be- 
tween same-momentum states with energy difference of 
hoj. This creates a circle in the FBZ where transitions be- 
tween the valence and conduction band are at resonance 
(see Fig. [IJ . We intend to use the modulation to reshuffle 
the spectrum such that the resulting valence band con- 
sists of two parts: the original valence band outside the 
resonance circle drawn in Fig. [T] and the original conduc- 
tion band inside the resonance circle, near the T point. 
From Fig. [I] we see that this indeed leads to the desired 
structure, with the reshuffled pseudospin configuration 
pointing south near the T point and north at large k- 
values (for M < 0). On the resonance circle, we expect 
an avoided crossing separating the reshuffled lower band 
from the upper band. 



FLOQUET TOPOLOGICAL INSULATOR IN A 
NON-EQUILIBRIUM (Cd, Hg)Te 
HETEROSTRUCTURE 



Let us next consider the Floquet problem in a 
Zincblende spectrum in detail. We add a time depen- 
dent field to the Hamiltonian ([5| 



V(t) = V • o-cos(wi), 



(9) 



where V is a three-dimensional vector, which has to be 
carefully chosen to obtain the desired result. It is con- 
venient to transform the bare Hamiltonian to a "rotat- 
ing frame of reference" such that the bottom band is 
shifted by fkj. This is achieved using the unitary trans- 
formation f7(k,i) = P + (k) + PL(k)e lLJt , where P±(k) = 
| / ± d(k) • <x are projectors on the of upper and lower 
bands of H(k). This results in the following Hamiltonian: 



Hj(Jt) = P+(k)e + (k) + P_(k) [e_(k) + u] 



UV{t)ll\ 
(10) 

where e± (k) are the energies corresponding to P± (k) . In 
the "rotating" picture, the two bands cross if w is larger 
then the gap M. The second term in the right-hand-side 



of Eq. ( 10 ) is the driving term, which directly couples the 



bands and has a time-independent component. 

The solution of Hj can also be given in terms of a 
spinor pointing along a unit vector, rik, which will play 
the same role as d(k) for the stationary iJ(k). fit will 
encode the topological properties of the FTI. 

Hi is solved by the eigenstates \ipf (k)), which for the 
values of momenta, k, away from the resonance ring are 
only weakly modified compared to the equilibrium, V — 
0, case. We define the vector n k = (ipj (k)\&\i/)J (k)), 
which characterizes the pseudospin configuration in the 
lower (— ) band of Hi. The vector nk is plotted in Fig. [2] 
for M/B < 0. Indeed, near the T point we see that 
nk points towards the south pole, and for larger values 
of k, the band consists of the original lower band, and 
therefore rik points towards the northern hemisphere for 
these k values. These two regimes are separated by the 
resonance ring, denoted by 7, for which uj — e + (k) — 
e~ (k) (the green curve in Fig. [l]) . 

The topological aspects of the reshuffled lower band 
depend crucially on the properties of fik on 7, which are, 
in turn, inherited from the geometric properties of the 
driving potential, encoded in V. These are best illus- 
trated by employing the rotating wave approximation, 
as we shall proceed to do below. An exact numerical 
solution will be presented in the next section. 

The driving field V(t) contains both counter-rotating 
and co-rotating terms. In the rotating wave approxima- 
tion it is given by 



W = P+ (k) ( V ■ &) P_ (k) + P_ (k) (V ■ &) P 4 



(k). 
(11) 




FIG. 2. Pseudospin configuration (blue arrows) and dis- 
persion of the lower band of Hi. Note the dip in the energy 
surface near k = 0, resulting from the reshuffling of the lower 
and upper bands of H(k). 



Next, we decompose the vector V as follows 

V= (v-d(k))d(k) + Vi(k). (12) 



A substitution in Eq. (Ill gives 



Vrwa = Vi ■ <t. 



On the curve 7 we have: 



Ak = -Vi(k)/|V_L(k)|. 



(13) 



(14) 



We can define a topological invariant C F similar to C 
in Eq. (JsJ) , by replacing d(k) with nk. In order for rik 
to have a non- vanishing C F , it needs to wrap around 
the unit sphere. A necessary condition is that on the 
curve 7, it forms a loop that winds around the north 
pole. Now, Vj_(k) lies on the plane defined by d(k) and 
V. For the values of k which are on-resonance, d(k) 
traces a closed loop on the unit sphere which encircles 
the north pole. Therefore, if the driving held vector V 
points to a point on the Bloch sphere which is encircled 
by this loop, Vj_(k) will also trace (a different) loop on 
the sphere which encircles the north pole, as illustrated 
by Fig. § 

Under the conditions stated above and with M < 0, 
the vector field starts from the south pole at T point 
and continues smoothly to the northern hemisphere for 
larger values of |k| while winding around the equator. For 
values of k further away from the curve 7, — d(k) as 

the driving field is off resonance there. The contribution 
of these k's to C F is therefore equal to their contribution 
to C. Therefore it is evident that C± — C± ± 1. Note 
that for M > 0, C F = C± T 1. 

A comment is in order regarding the time depen- 
dence of C F . Since the solutions to the time dependent 
Schrodinger equation are given by transforming back to 
the Schrodinger picture, 



d(k> 




FIG. 3. The geometrical condition for creating topological 
quasi-energy bands. The green arrow and circle depicts d(k) 
on the curve 7 in the FBZ for which the resonant condition 
holds. The red arrow and cone depicts Vj_(k) on 7. The 
blue arrow depicts the driving field vector V. As long as 
V points within the loop traced by d(k), the vector V_i_(k) 
winds around the north pole, which is indicated by the black 
arrow. 



the pseudospin configuration in the Brillouin zone of 
these solutions , 



n k (t) = (^-(M)|<x|<0-(k,i)), 



(16) 



will also depend on time. As long as Hi is non-degenerate 
in the FBZ, which implies that Ak(£) is well defined, C F 
will be quantized to an integer. As for C, also C F is a 
topological invariant which is robust to smooth changes 
in rik(t) which are not singular at any point in the FBZ. 
Therefore, C F does not depend on time, although the 
pseduospin configurations rik(i) do, and we can calculate 
using n k . 



c f 



NON-EQUILIBRIUM EDGE STATES 

One of the most striking results of the above consider- 
ations is the existence of helical edge states once the time 
dependent field is turned on. Below we demonstrate the 
formation of edge states in a tight binding model which 
contains the essential features of Eq. The Fourier 

transform of the spin-orbit coupling vector, d(k) in the 
corresponding lattice model is given by, cf., Eq. (TSl) , 



d(k) = (A 



sin k x , A sin k y , M — 4B 



- 2£?[cos k x + cos k y ]) . 

(17) 



\^(t,k)) = umf(k)), 



(15) 



We consider the above model with the time dependent 
field of the form Vq& z cos(wt) in the strip geometry, with 
periodic boundary condition in the x direction, and van- 
ishing boundary conditions at y = 0, L. 

We solve the Floquet equation numerically by moving 
to frequency space and truncating number of harmonics. 
The wave vector k x is therefore a good quantum number, 
and the solutions $(f) are characterized by e and k x . The 
quasi-energies for this geometry is displayed in Fig. [4] 
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(a) 



FIG. 4. Quasi-ene rgy spectrum of the Floquet equation |3| of 
the Hamiltonian 1 171, in the strip geometry: periodic bound- 



ary conditions in the x direction, and vanishing ones in the 
y direction. The driving field was taken to be in the z di- 
rection. The horizontal axis labels the momentum k x . The 
vertical axis labels the quasi-energies in units of \M\. Two 
linearly dispersing chiral edge modes traverse the gap in the 
quasi-energy spectrum. The parameters used are oj = 2.3jMj, 
|V| = A = \B\ = 0.2|M|. The inset shows the dispersion of 
the original Hamiltonian (171, for the same parameters. 



The quasi-energies of the bottom and top band represent 
modes which are extended spatially, while for each value 
of k x there are two modes which are localized in the y 
direction. 

As is evident from Fig. [4j the quasi-energies of these 
modes disperse linearly, e(k x ) oc k x , hence they are 
propagating with a fixed velocity. Consider a wave 
packet which is initially described by fo(k x ). From 
equation fcty we see that it will evolve into ip(t) = 
J dk x e i£< > ko =ftfo(k x )§l x (y,t), where denotes the quasi 
energy edge states with momentum k x . Clearly, this will 
give a velocity of (x) = J dk x \f(k x )\ 2 . 

In general, the solutions &e,k x (t) are time-dependent. 
An important finding is that the density edge modes are 
only very weakly dependent on time. This can be seen 
in Fig. [5j in which we plot the time dependence of the 
density profile of these modes. 



EXPERIMENTAL REALIZATION OF THE FTI 

To experimentally realize the proposed state, we 
need identify a proper time-dependent interaction in the 
HgTe/CdTe wells. Below we consider several options, 
of which the most promising one uses a periodic electric 




(b) 



FIG. 5. Density of edge mode as function of time, \<f>(y, t)\ 2 , 
(a) for k x = 0, and (b) for k x = 0.84, where the edge modes 
meet the bulk states. The horizontal axis display the distance 
from the edge, y, in units of the lattice constant, and the time 
in units of 2tt/uj. Only the density for the 20 lattice sites 
closest to the edge are shown for clarity. 

field, and the strong linear Stark effect that arises due to 
the unique band structure. 

Magnetic field realization - Perhaps the simplest real- 
ization of a time dependent perturbation of the form ^ 
is by a microwave-THz oscillating magnetic field, polar- 
ized in the z direction. The effect of Zeeman energies 
in thin Hg/CdTe quantum wells can be evaluated by re- 
calling that the effective model Q includes states with 
m.j = ±(1/2,3/2) in the upper and lower block respec- 
tively. This would result in an effective Zeeman splitting 
between the two states in each block [33]. The value 
for the g-factor for HgTe semiconductor quantum wells 
was measured to be g as 20 (23]. Therefore, a gap in 
the quasienergy spectrum on the order of 0.1K can be 
achieved using magnetic fields of lOmT. Bigger gaps may 
be achieved by using Se instead of Te, as its g-factor is 
roughly twice as large |25j . 

As can be seen by inspecting Eq. (12), the Chcrn 



numbers C F for each block in this realization depend only 
on the winding of the vector d(k). Therefore, the two 
blocks will exhibit opposite C F , resulting in two counter- 
propagating helical edge modes. As we explain in the 
next section, the counter-propagating edge modes cannot 
couple to open a gap in the quasi-energy spectrum, even 
though a magnetic field is odd under time reversal. 

Stress Modulation - Stress modulation of the quan- 
tum wells using piczo-clcctric materials, would lead to 
modulation of the parameter M in ([5]), and leads to two 
counter propagating edge states. 

Electric field realization - An in-plane electric field is 
perhaps the most promising route to the FTI, can pro- 
duce large gaps in the quasi-energy spectrum (compared 
the Zeeman case), and leads to robust co-propagating 
edge modes. The electric field is given by 



E = Re(£ • expiw£)iVk. 



(18) 



G 



Inserting this into Eq. ( 11 ), we get 



Vj_(k) = d(k)x(Re£-V k )d(k)-(Im£-V k )d(k). (19) 

As before, the vector field Vj_(k) is orthogonal to d(k), 
and again, we would like it to wind around the north 
pole. Now if we take £ = £(— ix + y) we get, expanding 
Eq.(19l to second order in k x ,k y , 



Vx(k) 



A(A 2 - 4BM)£ 



2 (^x 



(20) 



Evidently, the vector field Vj_ (k) winds twice around the 
equator. Therefore, for the above choice of £, the Chern 
numbers will be C± = ±2. We note that for the lower 
block will have C F — 0. Therefore, each edge of the 
system will have two co-propagating chiral modes, with 
spin in the sector of mj = 1/2,3/2. Naturally, a choice 
of £ = £(ix + y) will give C± = =f2 for the lower block 
and C F = for the upper block. The evolution with the 
oscillating electric field is topologically distinct from the 
trivial one, as the two co-propagating edge modes can- 
not become gapped. For HgTe/CdTe quantum wells with 
thickness of 58A51, we have |V(k)/£| ss 1mm , which 
leads to a gap in the quasi-energy spectrum on the or- 
der of 10K already for modest electric fields on the order 
of 1— which are experimentally accessible with powers 
< lmW. We note that decreasing the well thickness in- 
creases the value of |^4/M|, which can help achieve 
even larger gaps in the quasi energy spectrum. 



DISCUSSION 

In summary, we showed that the quasi-energy spec- 
trum of an otherwise ordinary band insulator irradi- 
ated by electromagnetic fields can exhibit non-trivial 
topological invariants and chiral edge modes. A real- 
ization of these ideas in Zincblcnde systems, such as 
HgTe/CdTe semiconducting quantum wells, can lead to 
Floquet Topological Insulators that support either co- or 
counter-propagating helical edge modes. The Floquet op- 
erators of these realizations belong, respectively, to sym- 
metry classes analogous to classes A (no symmetry) and 
All (time reversal symmetry with T 2 = —1) in [20] . 

The symmetry class of the Floquet Topological In- 
sulator indeed requires careful consideration when two 
counter-propagating edge states are present, as in the os- 
cillating magnetic-field realization suggested in the pre- 
vious section. In time independent systems, topological 
phases exhibiting counter-propagating edges are only dis- 
tinct from trivial phases under the restriction 1~HT~ 1 = 
H , where T is the anti-unitary time reversal operator sat- 
isfying T 2 = — 1. In the time-periodic case, the Hamil- 
tonian at any given time may not possess any symmetry 
under time reversal. Nevertheless, when the condition 



holds (for some r), the Floquet matrix of Eq. ^ satis- 
fies 7~&f (k)7~ — Hp(— k), where T is an anti-unitary 
operator which is related to T by T — V^TV, with V — 
<5k(-(r + r)/2), c.f. Eq. (||. Clearly, f 2 = -1. There- 
fore, under this condition, the quasi-energy spectrum 
consists of analogues to Kramer's doublets, which cannot 
be coupled by the Floquet matrix. The counter propagat- 
ing edge-modes are such a Kramer's pair, which, there- 
fore, cannot couple and open a gap (in the quasi-energy 



spectrum) under any perturbations satisfying Eq. (21 1 



We note that Eq. (21 ) holds for any Hamiltonian of the 
form T-L(t) = Ho + V cos (wt +(/>), with time reversal invari- 
ant Hq j and V having unique parity under time reversal, 
i.e., TVT^ 1 — ±V. An oscillating magnetic field, being 



odd under time reversal, therefore obeys Eq. (21) and 
leads to two counter propagating edge modes. 

Another important question that we did not touch 
upon is the non-equilibrium onset and steady states |27j 
of the driven systems. We emphasize that in the presence 
of the time dependent fields, response functions includ- 
ing Hall conductivity will be determined not only by the 
spectrum of the Floquet operator, but also by the dis- 
tribution of electrons on this spectrum. These in turn 
depend on the specific relaxation mechanisms present in 
the system, such as electron-phonon mechanisms [281 ES] 
and electron-electron interaction [30 . An interesting ob- 
servation here is that the average energies of the edge 
modes lie in the original bandgap. 

One way to minimize the unwanted non-equilibrium 
heating effects would be to use an adiabatic build-up of 
the Floquet topological insulator state, e.g. with the fre- 
quency of the modulation gradually increasing from zero. 
This should result, at least initially, in a filled lower Flo- 
quet band. Nevertheless, the scattering events will lead 
to a relaxation from the conduction band into the va- 
lence band of the original system and will always play 
role by producing mobile bulk quasi-particles. Perhaps 
this relaxation could be suppressed by restricting the cor- 
responding optical modes in the environment. An anal- 
ysis of the non-equilibrium states of the system will be 
the subject of future work. 
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TU{t)T- 1 =H(-t + T) 
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